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Neutrinoless double beta decay {[ifiOv) is tiie only realistic probe of the Majorana nature 
of the neutrino. In the standard picture, its rate is proportional to mee, the e — e element 
of the Majorana neutrino mass matrix in the flavor basis. I explore minimally allowed 
TOee values within the framework of mass matrix anarchy where neutrino parameters are 
defined statistically at low energies. Distributions of mixing angles are well defined by 
• the Haar integration measure, but masses are dependent on arbitrary weighting functions 

I and boundary conditions. I survey the integration measure parameter space and find that 

' for sufficiently convergent weightings, is constrained between (0.01 — 0.4) eV at 90% 

^ , confidence. Constraints from neutrino mixing data lower these bounds. Singular integration 

measures allow for arbitrarily small m^e values with the remaining elements ill-defined, but 
this condition constrains the fiavor structure of the model's ultraviolet completion. /3/30z^ 
bounds below mee ~ 5 x 10~^ eV should indicate symmetry in the lepton sector, new fight 
degrees of freedom or the Dirac nature of the neutrino. 
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^ , I. INTRODUCTION 

CN ' The observation of nonzero neutrino mass and mixing via flavor oscillations is the first terrestrial 

^ . evidence of physics beyond the Standard Model (SM) of particle physics. See [H, Q, [3| for a review of 

O ' neutrino physics. This discovery provides useful insight into the nature of new high scale phenomena 

^ I and introduces many important questions. Chief among these is the charge conjugation properties 

^ ■ (Dirac vs. Majorana) of the neutrino. First, one may enlarge the SM field content to include one or 

' more gauge singlets N (right handed neutrinos) and give the neutrino a Dirac mass similar to the 

00 , other fermions via coupling to the neutral SM Higgs boson: vNH^ . On the other hand, possessing 

no unbroken gauge quantum numbers, the neutrino may have a Majorana mass term which couples 
the neutrino to its charge conjugate v'^, and thus renders the neutrino equal to its own antiparticle. 
^ ■ This breaks all nontrivial global symmetries. In particular, a Majorana neutrino mass violates 

. lepton number by two units. The Dirac versus Majorana nature of the neutrino is an important 

issue that must be addressed experimentally. 

Neutrinoless double beta decay (/3/30i^) is an as yet unobserved Lepton Number Violating (LNV) 
process that would unambiguously identify the Majorana nature of the neutrino [1]. In fact, except 
in rare circumstances {3], is the only realistic hope of probing LNV in the near future 0, Q]- 

Given the importance of this process, it is useful to explore its theoretical expectations over a broad 
range of scenarios. In the standard picture, with no new light LNV degrees of freedom below the 
TeV scale, /3/30i/ proceeds primarily via Majorana neutrino exchange. In this case, its amplitude 
is proportional to mee, the e — e element of the neutrino mass matrix in the fiavor basis where 
the charged leptons are diagonal. While this relationship is generally spoiled by new physics it 
should still effectively hold for sufficiently small ruee values. While it is true that the smallest /3/30z/ 
rates may be dominated by high dimensional non-renormalizable operators there is an important 



feedback mechanism into rriee- This follows from the extended black box theorem ISH where a one to 
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one relationship is derived between rriee and the effective mll^ that governs PPOu in the vanishing 
limit, such that ruee = ^ mtl^ = 0. Thus, it is reasonable to assume that, when searching 
for small f3f30i' rates, it is enough to study the behavior of ruee- This is only loose motivation for 
the present analysis as it is impossible to extract the exact rate below which this reasoning holds 
without assuming properties of the neutrino mass's ultraviolet completion. In what follows I will 
assume the dominance of the light Major ana neutrino exchange mechanism parameterized by rUee- 
This is the most popular case. Additionally, since other new physics has yet to be discovered, 
it is the minimal mechanism currently implied by direct observation. Current experimental limits 
constrain the ffffOu half-life below ~ 10^^ years, corresponding to rUee < 0.35 eV at 90% confidence* 
Next generation experiments are poised to extend this reach by roughly an order 



of magnitude to niee < 0.05 eV [IJ, llfl, Il7l | 



While it is true that /jpOv rates are below current sensitivities^, it is possible that they will be 
discovered by the next round of experiments. In terms of the measured oscillation parameters, we 
are only beginning to explore the interesting range dominated by the atmospheric mass squared 
difference within the inverted and quasi-degenerate spectral hierarchies [l9|]. If /3/30i^ observation 
is "right around the corner," in which case TUee is relatively large, it is unlikely that its rate is 
suppressed by an approximate flavor symmetry. However, if bounds are pushed significantly lower, 
it is quite reasonable the a small niee is protected by an appropriate symmetry mechanism [20|]. It is 
natural to wonder how small it can be without the introduction of imposed mass matrix structure. 
To this end, it is instructive to take the "no structure" limit and consider niee bounds, assuming the 
anarchy hypothesis [2ll. l22l|. Similar reasoning was applied much earlier in the history of neutrino 
physics to probe the potential for larges mixing angles [23|]. In this scenario, the underlying neutrino 
mass model is sufficiently complicated, such that the mass matrix appears random at low energies 
in any basis. In other words, there is effectively no difference between the three light neutrino 



states. This leads to a distribution of observables that must be treated statistically. In [2]J, |2j], it 
was shown that the large mixing angles and small hierarchies of the neutrino sector are consistent 
with anarchy, provided ^13 is not too small, whereas the CKM matrix is inconsistent with anarchy, 
as expected by pure inspection of its structure-'- . In this analysis, the marginalized mixing angle 
distribution functions are well defined in terms of the Haar measure invariant under the U{3) 
group. It is not as straightforward to consider questions involving mass eigenvalues since one may 
include arbitrary U{3) invariant weighting functions into the integration measure [2l|| and boundary 
conditions. Here, I survey this added layer of ambiguity and derive the smallest allowed mass matrix 
element consistent with anarchy. 

In what follows, I analyze expectations for rfiee within the anarchy picture of neutrino mass 
generation. The goal is to determine how small/large one mass matrix element may be from the 
others within an anarchical framework. In section [ITl I introduce the formalism and notation em- 
ployed throughout the analysis. Using the Kolmogorov-Smirnov (KS) goodness of fit test, I scan 
the parameter space of measures defined by both simple polynomial and divergent U{3) invariant 
functions with "spherical" boundary conditions in Subsection III A[ Here, I also comment on modi- 



The translation bet-ween measured half-life and rriee is not straight for-ward, as it depends critically on isotope 
dependent nuclear matrix element calculations -where uncertainties currently range an order of mag nitude BHIlll- 
This is likely to improve -within the next several years. 

A positive signal -was reported by a subset of the HEIDELBERG-MOSCOW group -with a half-life near 1.19 x 10^^ 
years at 4.2(t confidence I neglect this observation in -what foUo-ws a-waiting conformation, except to point out 
that their extracted niee is -well accommodated by the anarchy model of neutrino mass. 

These claims -were questioned and explored by analysis discussed in [2^ and [2^ but the overall consistency bet-ween 
anarchy and current neutrino data still holds. 
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fications induced by the use of nontrivial boundary conditions. In Subsection III Bl I connect these 
results to the case of reahstic neutrino mixing. I conclude in Section IIIII with a summary of my 
results and comment on the impact of future experimental data. 

II. MASS MATRIX ANARCHY AND f3(30v 

I parameterize the complex, symmetric, three neutrino Majorana mass matrix as 

maf3 = maap = mrafse"*'"'^ , (11. 1) 

where the dimensionless complex parameters aa/3 = {a)ap define the structure of the matrix and 
are constructed to have a magnitude rap and phase (f)af3- The latter contains the familiar Majorana 
and Dirac phases of the neutrino mixing matrix in various linear combinations. Three of these 
six phases may be rotated away as unphysical with appropriate transformations, but are included 



in this analysis without loss of generality. See for example [27|, l28|] and references therein. A 
dimensionful factor of m is pulled out to carry the scale of neutrino masses. There is ambiguity 
in the factorization of m and r^^. To be concrete, I define m such that the resulting r^^ matrix 
elements have maximal magnitudes defined by boundary conditions subject to anarchy constraints. 
The average r^p values should be near unity. In other words, r^p is a generally 0{1) matrix up 
to some deviations described by the anarchy hypothesis. The overall neutrino mass scale m is 
inferred from experiment and included into the analysis by hand. Currently, m is bounded at 0.05 



eV from below by the atmospheric mass squared difference [l5|, l29(] and from above at roughly 1 eV 



by cosmological data [15||. In the remainder of this paper I refer to these as the hierarchical and 
quasi-degenerate limits, respectively. These are realized when the lightest mass eigenvalue squared 
is less than (hierarchical) or greater than (quasi-degenerate) the mass squared differences. In what 
follows, within the anarchy picture, one may only extract limits on the deviations of r^p from its 
average value, or equivalently, bound the ratio of mass matrix elements. 

Imposing the anarchy hypothesis implies that the matrix elements of Eq. (Ill.lj) are distributed 
randomly in any basis or, more precisely, the probability distribution of each Oq/3 is invariant under 
arbitrary unitary rotations. Starting from the diagonal mass basis, the flavor basis (as well as 
any other physical basis) is found by a random rotation, and one may calculate the probability 
distribution of any matrix structure. Notice that the obvious distinction of the diagonal mass basis 
renders it a very improbable structure which may be understood in terms of a sampled ensemble 
of matrices. The small chance of landing on it via random rotation does not preclude its existence. 

Assuming three light neutrinos, the condition of C/(3) invariance imposes strict conditions on 
the total distribution function G{aai3) universally associated with each element. This quantity may 
only be a function of 

det(a) = ei,-fcriir2,-r3fce^('^«+'^2^+'^3.) 

and the purely radial quantity 



tr(ata) = rj^ + r^^ + ^Is + 2?-?2 + 2r?3 + 2r, 



23- 



Working in these polar coordinates, the marginalized probability distribution of any single element's 
magnitude, say |aii| = rn = r, is obtained by integrating over the other radial coordinates and 
all phases, subject to some (also invariant) integration limits I = /[tr(a'^a), det(a)]. Thus, the 
r-coordinate distribution may be written in terms of the arbitrary functionals G and I, as 
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/■([tr(ata),dct(a)] ^ -, 

g(r) = / G tr(a"fa),det(a)J ?'Q/3C^'^a/3 JJ '^</'q/3- (11.2) 

At this point, G may be thought of as a weighting function in the integration measure. As a proba- 
bihty, g{r) is positive and normahzed as g{r)rdr = 1. This imphes that matrix elements do not 
have independent distribution functions, so that the probabihty of structures defined by magnitudes 
Tap cannot be expressed as a product oi g{rcip). Rather, one must calculate the generalized version 
Eq. (|II.2p . The cumulative probability distribution that the radial magnitude is less than some r 
is given by 

F{r) = I g{r')r'dr'. (II.3) 
Jo 

With this, the C confidence interval limits are found from 

max[F(r),l - F(r)] ^ = ^' ("•^) 

which yields two solutions interpreted as the extreme r values allowed at C confidence. One may 
trivially generalize this procedure to n neutrino states with no qualitative changes. It reflects the 
mechanism behind the KS goodness of fit test and yields the same results for one free variable. 



A. Simple Polynomial Measures 

Consider a one-dimensional marginalized probability distribution function, limited by / = 
Tr(a^a) < 3^. This choice is not unique but does make physical sense, as it implies that each matrix 
element lives within a three-dimensional sphere. This is clear upon rotation to the diagonal mass 
eigenbasis. I chose a non-unit radius for convenience to reflect the physical definition of the mass m, 
but this may be changed, provided a corresponding scaling of r. With this choice, the universal 0(1) 
matrix structures indicative of anarchy saturates the upper integration bound. This is easy to un- 
derstand, as sample volumes at large radii dominates that at smaller radii. Thus, r^/? G [0,3], with 
the most probable values naively expected near one. It is reasonable to assume the weighting func- 
tional may be expanded in a double Taylor series as G [tr(al^a), det(a)] = Cpgtr(a^a)^det(a)'^. 
Upon integrating over the complex phases it is clear that all nonzero powers of det(a) vanish by 
symmetry. However, this is only true with (pa/S independent integration limits. A nontrivial de- 
pendence on det(a) is possible within this framework but requires a departure from the physically 
motivated spherical boundary conditions. An example case of polynomial dependence reveals that 
such limits have little impact on the marginalized distribution function for small values of r, which 
is the primary concern of this analysis. 

Thus, it is enough to consider only linear combinations of 



Gp 



tr(ata), det(a)l = Tt{a^a)P = (r^ + r^^ + rfg + Irf^ + 2rl, + 2rl^Y . (II.5) 



In this case, it is simple to obtain the cumulative distribution function 
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FIG. 1: Anarchy allowed niee limits at 90% confidence as a function of p+6 defined by Eq. l|II.5|) (see text for 
details). Shaded region is total confidence interval while the blue (dash) and red (dash-dot) bounds indicate 
the hierarchical and quasi-degenerate neutrino mass spectra, respectively, niee as p ^ —6. Red and 
blue horizontal histograms show allowed upper and lower rriee bounds for a measure series expansion sample 
space for quasi-degenerate and hierarchical spectra, respectively. Current and future f3f30h' bounds are also 
shown. 



Fp{f) = (6+p)f^ 



+ ^ + 



5+p 4: + p 3+p 2+p 1+p {6+p){5+p){4+p){3+p){2+p){l+p) 

(11.6) 

via the procedure outlined in Eq. (|II.2p and Eq. (|II.3P in terms of the scaled magnitude f = r/3 
that lives in the unit interval. This is defined for all p > —6. It is easy to see that all other pole 
divergences cancel pairwise out of this expression. 

Combining Eq. (III.6|) and Eq. (|II.4p yields limits on the ratio of matrix elements r. Upon 
multiplication by the current experimentally allowed neutrino mass scale range 0.05eV < m < leV, 
one may obtain the two sided 90% confidence limits of rriee (or any other matrix element) as a 
function of p within the anarchy framework. This is shown in the gray contour of Figure [1] along 
with the current and future (3(30iy reach denoted by horizontal dotted and dashed lines, respectively. 
For convenience, the 90% allowed region is decomposed into the two extreme cases of hierarchical 
and quasi-degenerate mass spectra bounded by the blue dashed and red dotted-dashed curves, 
respectively. Given a particular neutrino mass scale, the allowed m-ee range spans scarcely an order 
of magnitude. Much of the breadth of the totally allowed grey region in Figure [T] comes from 
uncertainty in m. Therefore, an absolute neutrino mass measurement should significantly tighten 
the anarchy prediction of rriee. 

For p > —5 the curves of Figure [T] are relatively level yielding niee in the range (0.01-0.1) eV 
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and (0.05-0.5) eV for the hierarchical and quasi-degenerate spectra. This is easy to understand, 
as f < 1 and Eq. (|II.6P may be truncated at quadratic order^. Hence, Fp{r) — > ff^?^^ yielding 

a C confidence limit for r between 3y^(1 ifc C) iq^^^j^ - This approximation works very well for 
lower bounds but receives up to 5% corrections for the upper 90% confidence limits. At this 
point, it is instructive to consider the generalized n neutrino scenario for analogous weightings 
and boundary conditions I = Tr(a^a) < n^. In the small f = r/n approximation, the cumulative 
distribution function is given by Fp""^ w 2!!""^^^ { 2p+n"-i-n-2 } ' "^^ich yields confidence limits 

between (1 ± C) („_i)^(-^"!j.'^)(2p+n^+n) • large the matrix elements become independent with 
distributions bounded by ^Jl^C . 

I now consider general well behaved measure functions as series expansions of TT{a^aY weighted 
by arbitrary coefficients Cp, which by a trivial extension of Eq. (III.Gj) results in 

^ 00 

Here, one coefficient is always redundant, as it may be factored out of both the numerator and 
denominator and ultimately canceled. Due to this normalization condition, only those terms with 
comparable weightings having large destructive/constructive interference contribute to deviations 
from the limit contours of Figure [H This can be seen in the small f approximation F{r) « 

Z^^o / Sm=o 7^+6 • Uniformly scanning the parameters of Eq. (in.7p . truncated at p = 4 
within the range Cp S [—50,50], I sample the ruee bounds obtained from 10^'' sample functions. 
Results are shown in the blue and red horizontal histograms of Figure [U for the hierarchical and 
quasi-degenerate mass scales, respectively. These arbitrarily normalized distributions are highly 
dependent on the sampling procedure. Consequently, the upper and lower endpoints are the only 
useful quantities. Due to the small f approximate behavior, the distribution of upper bounds is 
much wider than that of the lower. The behavior of these histograms are in analogy with the work 
of [21| where random mass matrices were generated and studied using a linear measure function 
and cubic boundary conditions. In that case, mass eigenvalues were histogrammed as opposed to 
the combined quantity rriee- 

The behavior of Fp{r) for p — > —6 is easy to understand, but surprising from the mass matrix 
anarchy viewpoint. Approaching this lower limit, one finds that limp^_6 Fp{f) « f2(p+6) r^^i^^p^ 
6)lnf2, implying that the marginalized probability distribution g{r) — > d{r). This shows that at 
least one mass matrix element may be arbitrarily small, provided a sufficiently divergent integration 
measure. Furthermore, this must hold true in any basis obtained by a random rotation from the 
diagonal mass basis. With this, one finds the C confidence region bounded between e "(s+p) ^ 
which goes rapidly to zero, as seen in Figure [H It remains to check the behavior of the other, 
marginalized matrix elements in this limit. Given two magnitudes, r and s, one may calculate 
the marginalized probability distribution function gp{r,s) and the cumulative distribution function 
Fpir, s) parameterized by p in the integration measure of Eq. (|II.5p . In the limit p — > —6, when r 
is small, per the above argument, 

2 2 
Fp{r, s) « ^2(^^+6) -[p + af = Fp{r) - {p + 6)^. 



This is only true for p > — 5 as, in this limit, the quadratic terms of Eq. I|II.6P diverge and thus cancel. In this 
case one must truncate at quartic order. 
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Name 


Parameter Combination 


Value 


1(7 Uncertainty 


Am| 


9 9 

TOj — 


7.65 X 10-5 eV^ 


0.22 X 10-5 eV2 


Am^ 


\ml - mil 


2.40 X 10-3 eV^ 


0.12 X 10-3 


sin 9 s 


sin 9i2 


0.551 


0.017 


sin 6 a 


sin 6*23 


0.707 


0.046 


sin 6'_R 


sin 6*13 


0.1 


< 0.14 



TABLE I: Summary table of current neutrino results together with naming conventions and parameter 
definitions. Columns three and four list best fit central parameter values and la uncertainties. These were 
adapted from the global oscillation analysis of (s^l • The central values are used as input into the analysis of 
subsection III Bl 

When r is within its C confidence limits, defined by Fp{r) G (1 ± C)/2, then Fp{r, s) G (1 it C)/2 is 
also satisfied up to small perturbations. The coordinate s is virtually unconstrained. Thus, given 
p — > —6, at least one texture zero is guaranteed and all other matrix elements vary freely and 
are (almost) independent of the integration measure. Similar statements apply to other sufficiently 
divergent poles in the measure. For example, given G[Tr(a^a)] = (S^ — Tr(a^a))^ in the appropriate 
p limit, one radial magnitude goes to B while the rest are unconstrained by the theory. 



B. Mass Matrix Anarchy and Neutrino Mixing 



Thus far, I have explored the relation: Neutrino mixing data is consistent with the anarchy 
hypothesis... Given neutrino anarchy, what are the allowed mass matrix structures in an arbitrary 
basis? The allowed rriee range is then easily extracted. This is the proper treatment, but it does 
not address if/when the derived structures accommodate the neutrino data. Alternately, one might 
consider the anarchical distribution constrained by current oscillation phenomenology. This 
will not yield the same range as the previous analysis. As noted in [2j|, regarding the preferred 
value of the reactor mixing angle, given the consistency of neutrino data and anarchy, one would 
expect an expanded confidence region for unknown quantities such as the CP phases and 6*13. In 
other words, a large amount of (accidental) structure can appear among these parameters while 
still maintaining consistency with anarchy. In the standard parametrization, the composite mass 
matrix element rriee magnitude is given by 



m 



'12 cos + r2e 



^2 sin^ 



'13 



(11.8) 



where an overall mass scale m is factored out for consistency with the preceding analysis. Here 
and (pi make up the complex eigenvalues of the flavor basis mass matrix Eq. (Ill.ljl . diagonalized 
by the mixing angles 612, O23, 613 and 6 in the usual way. In this quantity, a single Major ana 
phase may be removed as unphysical so the resulting expression depends on only two phase linear 
combinations. The mixing parameters are distributed according to the Haar measure and the radial 
magnitudes according to the complex Majorana mass eigenvalue measure given in [2l||, weighted by 
an invariant function analogous to Eq. (jll.Sp . ruee is the summed convolution of these distributions. 
It is well known that, given current data, Eq. (Ill.Sj) can vanish, provided the normal neutrino mass 
ordering and particular relationships among the phases and mass scale [15|, l30|, [SlJ . One would then 
expect the allowed niee range to extend lower in the normal than the inverted hierarchy. 

I point out that the preceding analysis could have been done from this perspective, with sep- 
arated mass and mixing parameters, but was easier done in the flavor basis. Within a realistic 
neutrino mixing framework constrained by the data shown in Table [H the use of the convoluted 
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(III. 81) is convenient since we already know many best fit bounds from oscillation searches 
l29l |. Marginalizing the convoluted distribution over the unknown phases as well as the 
neutrino mixing angle uncertainties with a polynomial weighting function, I obtain the normalized 
iTiee cumulative distribution function. Thus, this quantity is the normalized term by term convolu- 
tion of angular distribution functions weighted by the Haar measure and the polynomial measured 
mass eigenstate distributions. The experimentally allowed mixing angle ranges may be substituted 
into this function. For simplicity, I consider the normal and inverted neutrino mass spectra sepa- 
rately in terms of only the lowest free mass eigenvalue with the others related to it by the measured 
mass squared differences. I find that for all non-singular weightings, niee > 4.4 x 10~^ eV and 
> 1-6 X 10~^ eV at 90% confidence for the normal and inverted hierarchies, respectively. For 
the inverted hierarchy, this bound coincides roughly with the smallest possible inverted ruee value 
obtained from Eq. (jll.SP evaluated at = 0. As expected, the normal hierarchy bound is well below 
those shown in Figured] due to cancellations induced by the marginalization over unknown phases. 
However, care must be taken when interpreting this result as one hierarchy may be preferred over 
another by anarchy. A numerical scan over an ensemble of mass eigenvalues, without imposing mass 
squared differences from oscillation data, reveals a general preference for the intermediate state to 
lie closer to the heavier than the lighter one. My results agree qualitatively with a similar scan 
done in [2l|]. Still, this slight effect does not suggest that anarchy favors the inverted hierarchy. 



which is defined in terms of mass squared differences. These results are relatively independent of 
the supplied weighting function, provided that it is sufficiently non-singular. 



III. CONCLUSIONS 

Within the framework of neutrino mass anarchy, the distribution of parameters and matrix ele- 
ments must be treated statistically. Unlike a study of neutrino mixing angles and phases, which de- 
pend only on the invariant Haar measure, an analysis of mass eigenvalues/matrix elements depends 
critically on arbitrary integration measures and boundary conditions. For well behaved measures, 
the value of any one matrix element may vary between about 0.01 eV and 0.4 eV at 90% confidence. 
This is well within the reach of future experiments and will be tightened with better knowledge 
of the overall neutrino mass scale. However, these bounds are expanded for sufficiently divergent 
measures, which may lead to a vanishing matrix element. This case also renders the remaining 
matrix elements undefined, and thus removes almost all predictability from the theory. I conclude 
that arbitrarily small matrix elements are allowed within the anarchy framework, but this seems 
counter-intuitive, given the universal mass matrix structure assumed in the original formulation of 
the theory. 

This is more puzzling when put in the framework of realistic neutrino data, niee = implies 
the normal hierarchy and practically forces a — r like symmetry among the remaining elements 
[^,[33, 34|. While still consistent with anarchy (the divergent weighting imposes the texture zero and 



frees all other elements to vary almost unconstrained to yield the correct mass spectra and parameter 
relationships), this seems like a lot of accidental structure for a structureless matrix! The resolution 
lies in the ultraviolet completion of the theory that produces the integration measure. That is, 
the underlying theory yields well defined mass parameters, but by our limited knowledge of its 
complicated structure, this information is communicated to low energies as an ensemble of possible 
choices that must be treated statistically. The nature of possible model classes select the weighting 
function for us. An anarchy guaranteed texture zero implies a mechanism selecting those textures 
from the ensemble of ultraviolet completions. This is inconsistent with the anarchy principle from 
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construction, as it would require a corresponding symmetry mechanism and/or parameter fine 
tuning. Of course, without knowing the details of neutrino mass generation it is impossible to 
gauge or select uncomfortable levels of tuning or when flavor structures must arise. Still, one would 
not expect the metric power law parameter to venture too close the the —6 singularity. To get a 
handle on this, it is reasonable to assume a 10% or greater deviation which implies rriee > 5 x 10"^ 
eV at 90% confidence. [ifiQv bounds below this level should indicate nontrivial structure in the 
lepton flavor sector, new light LNV degrees of freedom 
neutrino. 



35, 



36, 



37, 



38j or the Dirac nature of the 
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